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Abstract 

Salmonella enterica subsp. enterica serovar Enteritidis (5. Enteritidis) was responsible for a worldwide pandemic 
during the 1980s and 1990s; however, changes in the dominant lineage before and after this event remain 
unknown. This study determined 5. Enteritidis lineages before and after this pandemic event in Japan using 
multilocus sequence typing (MLST). Thirty 5. Enteritidis strains were collected in Japan between 1973 and 2004, 
consisting of 27 human strains from individual episodes, a bovine strain, a liquid egg strain and an eggshell strain. 
Strains showed nine phage types and 17 pulsed-field profiles with pulsed-field gel electrophoresis. All strains had 
homologous type 1 1 sequences without any nucleotide differences in seven housekeeping genes. These MLST 
results suggest that 5. Enteritidis with the diversities revealed by phage typing and pulsed-field profiling has a 
highly clonal population. Although type 11 5. Enteritidis may exhibit both pleiotropic surface structure and pulsed- 
field type variation, it is likely to be a stable lineage derived from an ancestor before the 1980s and/or 1990s 
pandemic in Japan. 



Findings 

Salmonella enterica subsp. enterica serovar Enteritidis 
(«S. Enteritidis) was responsible for a worldwide pandemic 
during the 1980s and 1990s; however, changes in the 
dominant lineage before and after this event remain 
unknown. It is also difficult to explain why multiple clones 
of S. Enteritidis (defined by phage types, (PTs)) infecting 
chicken reproductive tissues emerged simultaneously in 
geographically separate countries during this pandemic. 
For example, PT 4 emerged in some European countries 
during the late 1980s, including England and Wales [1], 
France [2] and Germany [3], while PT 6 dominated Den- 
mark [4]. PT 8, PT 13a and PT 13 were most common in 
the USA [5] and in Japan, PT 1 and PT 4 were dominant 
[6]. To understand the simultaneous emergence of this 
pathogen, it is necessary to deduce whether these multiple 
lineages co-instantaneously acquired the ability to colonize 
chicken reproductive tissues or a single clone that 
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possessed this ability evolved into multiple lineages before 
the pandemic. 

Multi-locus sequence typing (MLST) can be used to 
accurately identify bacterial lineages [7,8]. The genetic dis- 
tance between two strains can be quantitatively estimated 
as allelic differences in the nucleotide sequences of house- 
keeping or virulent genes among bacterial strains [8]. 
MLST data have been used in evolutionary and population 
analyses that estimate recombination and mutation rates 
and also in the investigation of evolutionary relationships 
among bacteria classified within the same genus [8]. Using 
MLST, our study aimed to determine the lineages of 
S. Enteritidis before and after the pandemic during the 
1980s and/or 1990s in Fukuoka Prefecture, Japan. 

Salmonella strains used in this study are listed in 
Table 1. Human strains from individual food-borne 
disease outbreaks in 1973-1981 were isolated in the 
Kanagawa Prefectural Institute of Public Health. Other 
human strains (isolated from individual outbreaks), the 
liquid egg strain and the eggshell strain were isolated 
from different samples at the Fukuoka Institute of 
Health and Environmental Sciences. The cattle strain 
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Table 1 Phage type, pulsed-field gel electrophoresis profile and sequence type of Salmonella enterica subsp 

Strain No. Isolation Date Origin Foods that caused disease Phage Type PFGE Profile ST a 
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enterica serovar Enteritidis strains in Japan from 1973 to 2004. 
a ST, strain type. 
b NT, not tested. 



was provided by the Livestock Hygiene Service Center. 
All strains were isolated in Japan between 1973 and 
2004 (Table 1). We selected these strains to assess long- 
term alternation or rotation of S. Enteritidis genotypes. 
The study has been approved by the Fukuoka Institute 
of Health and Environmental Silences Ethics Committee. 

Housekeeping genes, thrA (aspartokinase+homoserine 
dehydrogenase), purE (phosphoribosylaminoimidazole car- 
boxylase), sucA (alpha ketoglutarate dehydrogenase), hisD 
(histidinol dehydrogenase), aroC (chorismate synthase), 
hemD (uroporphyrinogen III cosynthase) and dnaN (DNA 
polymerase III beta subunit) were used for MLST charac- 
terization. Each overnight-culture of the strains grown in 
LB broth (Difco Laboratories, Sparks, MD., USA) was cen- 
trifuged at 12,000 rpm for 5 min before the supernatant 
was discarded. The pellet was resuspended in normal sal- 
ine (1 ml) and centrifuged at 12,000 rpm for 5 min. The 



supernatant was discarded and the pellet was resuspended 
in 200 ul of 0.1 x TE buffer before being boiled for 10 min 
and centrifuged at 12,000 rpm for 10 min. The superna- 
tant was used as a template for the primary PCR. Primers 
for the primary PCR were as follows: thrA F 5'-GTC ACG 
GTG ATC GAT CCG GT-3'; thrA R 5'-CAC GAT ATT 
GAT ATT AGC CCG-3'; purE F 5'-ATG TCT TCC CGC 
AAT AAT CC-3'; purE R 5'-TCA TAG CGT CCC CCG 
CGG ATC-3'; sucA F 5'-AGC ACC GAA GAG AAA 
CGC TG-3'; sucA R 5'-GGT TGT TGA TAA CGA TAC 
GTA C-3'; hisD F 5'-GAA ACG TTC CAT TCC GCG 
CAG AC-3'; hisD R 5'-CTG AAC GGT CAT CCG TTT 
CTG-3'; aroC F 5'-CCT GGC ACC TCG CGC TAT AC- 
S'; aroC R 5'-CCA CAC ACG GAT CGT GGC G-3'; 
hemD F 5'-ATG AGT ATT CTG ATC ACC CG-3'; hemD 
Fl 5'-GAA GCG TTA GTG AGC CGT CTG CG-3'; 
dnaN F 5'-ATG AAA TTT ACC GTT GAA CGT GA-3'; 
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dnaN R 5'-AAT TTC TCA TTC GAG AGG ATT GC-3\ 
Amplifications for all genes were carried out with 2 \A 
DNA template, 2 \A of 2.5 mM (each) deoxynucleoside tri- 
phosphates, 5 ul of 5 x PrimeSTAR Buffer (5 mM Mg 2+ ) 
(Takara Bio Inc., Otsu, Japan), 6.25 pmol of primer, and 
0.625 U of PrimeSTAR HS DNA Polymerase (Takara Bio 
Inc.) in 25 \i\ reaction mixtures. Reaction conditions 
involved 30 cycles of: (A) denaturation at 98°C for 10 s, 
(B) primer annealing at 55°C for 10 s, and (C) extension at 
72°C for 1 min. Reactions were performed using a PCR 
Thermal Cycler SP TP 400 (Takara Shuzo, Co., Ltd., 
Kyoto, Japan). PCR products were purified using Microcon 
(Millipore Corporation, Bedford, Mass., USA) and 
sequenced using an ABI Prism BigDye Terminator v3.1 
Cycle Sequencing Kit and an ABI 3730 Genetic Analyzer 
(Applied Biosystems Japan Ltd., Tokyo, Japan). Primers for 
sequence reactions were as follows: thrA-Y 5'-ATC CCG 
GCC GAT CAC ATG AT-3; thrA-R 5'-CTC CAG CAG 
CCC CTC TTT CAG-3'; purE-Y 5'-CGC ATT ATT CCG 
GCG CGT GT-3'; purE-R 5'-CGC GGA TCG GGA TTT 
TCC AG-3'; sucA-Y 5'-AGC ACC GAA GAG AAA CGC 
TG-3'; sucA-R 5'-GGT TGT TGA TAA CGA TAC GTA 
C-3'; hidD -F 5'-GTC GGT CTG TAT ATT CCC GG-3'; 
hisD-R 5'-GGT AAT CGC ATC CAC CAA ATC-3'; 
aroC-Y 5'-GGC ACC AGT ATT GGC CTG CT-3'; aroC- 
R 5'-CAT ATG CGC CAC AAT GTG TTG-3'; hemD-F5'- 
GTG GCC TGG AGT TTT CCA CT-3'; hemD-R5'-GAC 
CAA TAG CCG ACA GCG TAG-3'; dnaN-Y 5'-CCG 
ATT CTC GGT AAC CTG CT-3'; dnaD-R 5'-CCA TCC 
ACC AGC TTC GAG GT-3\ Sequences were edited and 
the complementary sense and antisense fragments were 
aligned using the SeqManll program in the Lasergene 
software package (DNASTAR, Madison, WI, USA). 

Sequences were submitted to the MLST database website 
(http://mlst.ucc.ie/) and assigned existing or novel allele 
type numbers and sequence type numbers defined by the 
database. Composite sequence types (STs) were assigned 
based on the set of allele types derived from each of the 
seven loci. Allele sequences for each strain were then con- 
catenated in the order aroC - dnaN - hemD - hisD - purE 
- sucA - thrA for a final composite length of 3,336 bp. Con- 
catenated sequences were aligned using ClustalW software 
[9] and a phylogenetic tree was constructed using the 
neighbor-joining method (Center for Information Biology 
and DNA Data Bank of Japan) to compare to other pub- 
lished STs. Phylogenetic analyses were performed using 
NJplot (http://pbil.univ-lyonl.fr/software/njplot.html). 

Phage typing was conducted on S. Enteritidis strains 
recovered from samples in this study according to the 
method of Ward et al. [10] at the National Institute of 
Infectious Disease. 

Pulsed-field gel electrophoresis (PFGE) was conducted 
as described in our previous report [11]. Strains were 
considered different if variation in one or more DNA 



bands was observed. Using FPQuest Software (Bio-Rad 
Laboratories Inc. Hercules, CA, USA), similarity and 
cluster analyses were performed with Dice coefficients 
of similarity and an unweighted pair group method 
using average linkage, respectively. 

Phage typing identified nine PTs (Table 1). PT 1 was 
most dominant, while PT 3, PT 7 and PT 36 were indivi- 
dually detected in only one strain for each. For PFGE ana- 
lysis, 30 strains generated 17 pulsed-field profiles (PFPs) 
after restriction with Blnl (Table 1 Figure 1). PFP 1 was 
most dominant (11 strains), followed by PFP 10 (three 
strains). All other PFPs were only found in a single strain. 
S. Enteritidis with PFP 1 consisted of PT 1, PT 4, PT 14b 
and PT 36 strains. 

MLST analysis did not detect any nucleotide differ- 
ences among the 30 S. Enteritidis strains for seven genes 
(thrA: 501 bp, pure: 399 bp, sucA: 501 bp, hisD: 501 bp, 
aroC: 501 bp, hemD: 432 bp, and dnaN: 501 bp). All 
strains were characterized as ST 11 according to the 
MLST database (Table 1). The sequences of ST 11 are 
available in the MLST database http://mlst.ucc.ie/. ST 
11 consisted of strains with variable PTs and PFPs. 
According to concatenated sequences analysis, ST 11 
belonged to a clade that consisted of ST 11, ST 136, 
and ST 183 (Figure 2). 

For the first reported time, our study characterized all 
tested S. Enteritidis strains between 1973 and 2004 as a 
single ST (ST 11) and several PTs and PFPs. Two main 
interpretations of these results include: (A) Rabsch and 
colleagues speculated that S. Enteritidis acquired a niche 
in chicken reproductive tissues before the 1980s [12]; 
(B) disease-associated S. Enteritidis had a higher clonal 
population in Fukuoka Prefecture, Japan. These hypoth- 
eses are further discussed below. 

Persistence of S. Enteritidis clonality in Japan supports 
the idea that the ST 11 lineage had possibly acquired a 
niche within chicken reproductive organs before the pan- 
demics (1980s and 1990s). It is reported that the process 
of chicken oviduct colonization of S. Enteritidis is complex 
and depends on many factors including fimbriae, flagellae, 
lipopolysaccharide, cell wall structure and stress tolerance 
[13]. Until now, it has not been clear as to why S. Enteriti- 
dis is specifically implicated in egg contamination [14]. 
Housekeeping genes potentially play a less important roll 
in chicken oviduct colonization than genes of the above- 
mentioned factors. However, the absence of changes in 
the housekeeping genes suggests that the ST11 lineage 
survived without any substantial genetic changes being 
produced by the fundamental processes of mutation, 
recombination, selection, drift, and migration [15]. Insuffi- 
cient information was available for strains isolated before 
the pandemic because only three isolates were tested. 
Although ST 11 S. Enteritidis may produce some strain 
variation in pleiotropic surface structures [16] defined by 
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Figure 1 Pulsed-field gel electrophoresis analysis of fingerprints from genetically defined Salmonella enterica subsp. enterica serovar 
Enteritidis strains. Numbers on the right of the dendrogram indicate pulsed-field profiles (PFPs) and the phage types (PTs) associated with each 
PFP. The scale indicates the percentage similarity, as determined using Dice coefficients. PFPs 1 and 10 include some strains. The center 
diagrams are the PFPs of Salmonella enterica subsp. enterica serovar Enteritidis strains with Blnl digestion. 



PTs (such as PT 1, 3, 4, 6a, 7, 9, 14b, and 36) and geno- 
types defined by PFPs, it is likely to be a stable lineage that 
was derived from an ancestor before the pandemic in 
Japan. Our study therefore potentially supports the specu- 
lative assumptions made by Rabsch et al. [12]. 

Our study also suggested that food-borne disease- 
associated S. Enteritidis has a relatively higher clonal 
population than other Salmonella. In another study, 81 
isolates of S. Newport from humans, food animals, and 
retail foods showed 12 STs [17], while in our present 
study, S. Enteritidis populations showed a single ST line- 
age. As all typing systems are influenced by the tested 
population [18], we can only discuss human and food 
related strains in Fukuoka Prefecture, Japan, rather than 
worldwide S. Enteritidis population genetics. However, 
according to the international MLST database, among 
the several distinct STs of S. Enteritidis in the world, 
only ST 11 consists of human strains. These results 



imply that human and food related-S. Enteritidis may be 
a higher clonal population than S. Newport. 

MLST is a useful method for analyzing the core genes 
of pathogens that are of public health importance. 
MLST can be used to analyze conserved core genes, 
like housekeeping genes, which are universally con- 
served and provide fundamental genetic information 
[19]. Our study has revealed that MLST is suitable for 
evolutionary, basic relatedness or population studies, 
however, in contrast with PFGE, this method is not 
effective for epidemiological purposes. This is because 
MLST cannot analyze a diverse range of Salmonella 
strain specific genetic information that may be asso- 
ciated with niche adaptation. However, genetic variation 
in pathogen populations is a major barrier to disease 
control [20], and therefore, it is pertinent that we con- 
tinue to survey the S. Enteritidis lineage in Japan using 
MLST. 
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Figure 2 Phylogenetic tree of Salmonella enterica subsp. enterico serovar Enteritidis and Salmonella enterica subsp. enterica serovar Typhi 
according to STs from the MLST database (http://mlst.ucc.ie/mlst/dbs/Senterica), showing the difference between ST 1 1 and others. All strains 
from our study were ST 1 1. Allele sequences for each strain were concatenated in the order aroC - dnaN - hemD - hisD - pur\L - sucA - thrA for a 
final composite length of 3,336 bp. Sequences were aligned using the ClustalW method [9] and phylogenetic analyses were performed using 
NJplot through the Web Search and Analysis Service System provided by DNA Data Bank of Japan (Center for Information Biology and DNA 
Data Bank of Japan). All strains were designated in ST1 1. The scale represents nucleotide substitutions per site. 
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